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We study the physics of the superconducting variant of Weyl semimetals, which may be realized in 
multilayer structures comprising topological insulators and superconductors. We show how super- 
conductivity splits each Weyl node into two. The resulting Bogoliubov Weyl nodes can be pairwise 
independently controlled, allowing to access a set of phases characterized by different numbers of 
bulk Bogoliubov Weyl nodes and chiral Majorana surface modes. We analyze the physics of vortices 
in such systems, which trap zero energy Majorana modes only under certain conditions. We finally 
comment on possible experimental probes, thereby also exploiting the similarities between Weyl 
superconductors and 2-dimensional p + ip superconductors. 

PACS numbers: 



I. INTRODUCTION 

The discovery of topological insulators has stimulated 
a broad inquiry into topological features of electronic en- 
ergy bands. Such features are present not only in fully 
gapped systems but also in gapless ones. Of particu- 
lar recent interest are semimetals (zero gap semiconduc- 
tors) with Fermi points, where conduction and valence 
bands touch at isolated momenta in the Brillouin zone. 
In three dimensions, a linear touching between two non- 
degenerate bands is a Weyl point and is completely 
robust to all perturbations which do not break trans- 
lational symmetry. Such Weyl nodes are predicted to 
lead to a variety of measurable consequences, includ- 
ing unusual surface states whose Fermi surface is open 
("Fermi arcs"), unusual Hall effects, and other unusual 
transport features, and have been studied in a number 
of systems Generically, Weyl points require a system 
with strong spin-orbit coupling, and in addition the con- 
dition that the bands be non-degenerate requires that 
at least either inversion or time-reversal symmetry be 
broken.—^ This might occur in a bulk material through 
magnetic order or a non-centrosymmctric crystal struc- 
ture, but Weyl points can also be engineered. In particu- 
lar, an appropriate superlattice of alternating normal and 
topological insulators has been shown to display Weyl 
points£ii2 

Weyl fcrmions have been discussed extensively in the 
context of the A phase of 3 He, which also exhibits Weyl 
quasiparticles. 3 He also supports the B phase, in which 
quasiparticles are fully gapped, but which nevertheless 
possesses interesting topological properties^ In this pa- 
per, we explore the connection of Weyl semimetals to 
Weyl and topological superconductors, and more gener- 
ally the effects of superconductivity on Weyl semimet- 
als. Specifically we consider another class of engineered 
structures, in which the normal insulator of the afore- 
mentioned superlattice is replaced by a (conventional, 
s-wave) superconductor. Simple arguments show that 
only the Weyl semimetal produced by time-reversal sym- 
metry breaking, and not the one produced by non- 



ccntrosymmetry, leads to non-trivial superconducting 
states. Focusing on the former, we find that a vari- 
ety of superconducting phases, with and without Weyl 
points, and with varying topological features, can be 
tuned, depending upon the degree of time-reversal sym- 
metry breaking and the magnitude of the superconduct- 
ing proximity effect upon the topological insulating layers 
of the superlattice. 

Each of these phases may be characterized in a number 
of ways, which we discuss in the main text. In the bulk, 
they may be parametrized by the number, location, and 
chirality of Weyl points. At surfaces, depending upon 
both the phase and the surface orientation, one or several 
branches of chiral Majorana states may be present, and 
extend over a varying range of momenta in the surface 
Brillouin zone. Gapless Majorana states may also be 
present at vortex cores, again depending in detail upon 
vortex orientation and phase. These Majorana modes are 
relatives of those proposed for use in quantum computing 
in two dimensional topological superconductors. 

The paper is organized as follows. In Sec. [Til we dis- 
cuss how Weyl superconductors can be engineered in su- 
perlattices of superconductors and topological insulators, 
and derive the corresponding model. Sec. Illll is then de- 
voted to the characterization of this Hamiltonian. We 
analyze how superconductivity acts on Weyl electrons 
in the bulk and discuss the related topological surface 
physics. Based thereon, a topological phase diagram is 
constructed. In Sec. IIV1 we analyze the physics of vor- 
tices, which under certain conditions can bind Majorana 
zero modes. We close with some proposals for experimen- 
tal signatures based on thermal and electrical transport, 
which arc given in Sec. |Vl 



II. SUPERLATTICE 

A Weyl semimetal can be understood as an intermedi- 
ate phase between a normal insulator (NI) and a topolog- 
ical insulator (TI), arising due to a perturbation of the 
transition between the twoJ^ One pathway to engineer 
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FIG. 1: The proposed heterostructure for an experimental re- 
alization of Weyl superconductors. Magnetically doped layers 
of a topological insulator (TI) are alternated with layers of an 
s-wave superconductor (SC). The period of the superlattice is 
d. The arrows in the TI layers depict their magnetization, 
which is along the superlattice axis. 



Weyl semimetals is thus the stacking of layers of topo- 
logical and normal insulators^ In the very same spirit, 
a Weyl superconductor arises upon alternating thin lay- 
ers of topological insulators and standard s-wave BCS 
superconductors (SC). In such a structure, sketched in 
Fig. [TJ the proximity effect induces superconductivity in 
the surface states of the TI layers. 

As remarked in the Introduction, to realize a Weyl 
scmimetal requires breaking of cither time reversal or 
inversion symmetry. In the bulk of this paper we fo- 
cus on the time reversal symmetry breaking case, as it 
leads to much more non-trivial results in the presence 
of superconductivity. Indeed, in Appendix [AJ we show 
that when inversion symmetry is instead broken, while 
time reversal is preserved, superconducting proximity ef- 
fect leads directly to a gapped, trivial phase. Specifi- 
cally we model the topological insulator layers by a single 
Dirac node per surface, with an imposed Zccman splitting 
which may be considered to arise from an exchange cou- 
pling to randomly distributed magnetic impurities fer- 
romagnetically polarized perpendicular to the TI layers, 
while the coupling between the TI surface modes and the 
magnetic field of the impurities is neglected as usual in 
such systems i^2r— We assume the individual layers to be 
thin, such that sample is globally phase coherent. Fur- 
thermore, neighboring surface layers are tunnel coupled. 
Longer range tunneling is assumed to be negligible. 



A. Model Hamiltonian 

Working in units of h = 1 , wc model the system by the 
Hamiltonian 



H = E 4 H ^ c t + H sc 

t-^i k±i J k±j 



(1) 



Hij = Vjp t z (z x a) ■ k± Sij +mcr z 5ij (2) 



h sc =Ya c i op V°p t + c bot t c bot. t + hx _ 

k± ,i 

(3) 



where eg 
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V fcj_Ti' k ± li' fexti' k ± li J 



comprises an- 
nihilation operators for electrons of spin up and down 
in the top and bottom surfaces of layer i with in-plane 
momentum kj_ . The unit vector along the perpendicular 
axis is z. The Fermi velocity of the Dirac nodes is vf, 
for simplicity considered to be the same on each surface, 
and Pauli matrices a act on the real spin. The addi- 
tional pseudo spin for the top/bottom surface degree of 
freedom denoted by the Pauli matrices f. The Zeeman 
mass (half the Zeeman splitting) is m, the tunneling be- 
tween top and bottom surface of the same TI layer is 
denoted by ts, and the tunneling between different TI 
layers is tjj . The proximity induced superconductivity is 
characterized by A = |A|e lv , with (p being the globally 
coherent superconducting phase. 

Wc proceed by Fourier transforming the Hamiltonian 
along i, where the superlattice constant is d. After a 
canonical transformation 



a ± -> r z a ± 



r ± -> cr z r ± , 



(4) 



and subsequent diagonalization in the f subspace, the 
Hamiltonian reads 



SC,l 



H± =v F (z x a) ■k + M±(k z )a z , 



(5) 



(6) 



M±(k z ) =m± Jt 2 s + t 2 D + 2t s t D cos (k z d) , (7) 



(8) 



where cg ± 



(c- , c r ) T is now composed of the appro- 



priate eigcnoperators resulting from the diagonalization 
in the r-subspace, and k being the 3-dimensional momen- 
tum. Because the proximity induced superconductivity 
does not mix the two f sectors, we can analyze the cor- 
responding subspaccs separately. 



B. Normal state 

A Weyl node corresponds to the vanishing of an 
eigenenergy of ([S]) for one momentum. In the case with- 
out superconductivity analyzed in Ref. [H, and assuming 
without loss of generality that m > and is/to > 0, two 
Weyl nodes of opposite chirality can appear in the spec- 
trum of H-. They are located at k 
with 



(0,0,7r/d±/c ) 



ko = — arccos 
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as long as the condition 
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(9) 
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is fulfilled. For to 2 < (tg — to) 2 , H- describes a trivial 
insulator, while to 2 > (tg + to) corresponds to a quan- 
tum anomalous Hall insulator. The Hamiltonian W+, on 
the other hand, always describes a trivial insulator. It is 
adiabatically connected to the case m = that is topo- 
logically trivial. 



For m > |A|, this subspace is adiabatically connected 
to the topologically trivial case to = |A| = 0, and Weyl 
nodes can only appear in the t z = — I sector. By anal- 
ogy to the normal case, we find that for m > |A|, the 
spectrum of Eq. (fl"2)) has up to 4 Bogoliubov Weyl nodes 

of pairwise opposite chiralities at k = (0,Q 1 ir/d± k±) T 
with 



C. Superconducting state 

For the superconducting case, wc keep to > and 
tg/tD > 0, although a different choice does not change 
our results qualitatively. We start by analyzing the sub- 
space corresponding to r z = —1, which in the normal 
case potentially exhibits Weyl nodes. Technically, su- 
perconductivity is taken into account by introducing a 
particle-hole pscudospin on which the Pauli matrices k 
act, as well as the corresponding Nambu spinors. Using 
^ = (%_,c^_,cl u _,cl^_) T , the r z = -1 sector of 
the Hamiltonian can be recast into the form 



v F (z x a) ■ k + M-(k z )a z ) tg (11) 
. 1 
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Diagonalization of the k subspace yields 



H- = - V $t W 

9 Z-^i k.n 



where 

H ±A 
M ±A 



-v F (zxa)- k + M^(k z )a 



(12) 



(13) 



(k z ) =(m± |A|) - \lt% + t 2 D +2 tgt D cos (k z d) , 

(14) 
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In the basis of Bogoliubov quasiparticles and /g, the 
Hamiltonian thus takes a similar form to a normal 
Weyl semimctal upon replacing m — > m ± |A|. 

The subspace corresponding to t z = +1 can be ana- 
lyzed in the same way, which leads to a Hamiltonian 

H% A =v F (z x a) ■ k + AI± A (k z )a z , (16) 

Mf A (k z ) = (m ± |A|) + yJt 2 s + t 2 D +2t s t D cos(M) . 

(17) 
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(m±|A|) 2 -(t s -* D ) ; 
if the respective conditions 



m c i < m ± | A | < to 



c2 



(18) 



(19) 



are fulfilled. For m ± |A| < m cl , the respective mode 
is adiabatically connected to the case to = |A| = and 
thus topologically trivial. 

If to < |A|, each f sectors contains one topologically 
trivial mode as well as one mode that potentially has 
Weyl nodes. The latter now exist in the range 



TOcl < |A| ± to < to c2 , 
at the same momenta k = (0, 0,7r/d± fc±) T with 



(20) 



and the topologically trivial regime is corresponds to 
| A| ± to < m c \. 



III. CHARACTERIZATION OF THE 
ACCESSIBLE PHASES 

Having recast the Weyl superconductor Hamiltonian 
into a more convenient form, we will now analyze it in 
detail. For simplicity, we focus on the case to > |A|, 
when all of the interesting physics happens in H- defined 
in Eq. (fT2"j) . As the discussion in the last section implies, 
the case to < |A| follows upon interchanging the roles of 
to and |A| and considering different subspaces of the full 
Hamiltonian. 

We first recall some results for the limiting case | A | — > 
0, in which a normal Weyl semimetal is recovered (see 
Fig. [2ta)).— If the Zeeman mass to is small, to < m c i, 
the system is a topologically trivial insulator. When m 
reaches m c i , two Weyl nodes of opposite chirality appear 
at K = (0, 0, Tr/d) T . Close to these points, k = K + q, 
the dispersion is roughly given by 



E rj ±vp a ■ q 



(22) 



where ± defines the chirality of the node. 

Upon increasing to, the Weyl nodes move in oppo- 
site directions along the fc z -axis and to the momenta 
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FIG. 2: Evolution of the masses MZ A (upper curve) and 
AlZ A (lower curve) defined in Eq. (| 14[) upon increasing |A|. 
For |A| = and m c i < m < m C 2, the system has two Weyl 
nodes of chiral electrons, located at the sign changes of M_ A . 
With superconductivity, each Weyl nodes splits into two Bo- 
goliubov Weyl nodes of equal chirality and opposite particle- 
hole symmetry. Their separation grows with increasing |A| 
from subfigures (a) to (d). 



k = (0, 0, n/d ± k A ~°) . For fixed k z , the combined 
Hamiltonians H^ A ~° describe a gapped 2-dimensional 
Dirac electron. The mass of the latter changes sign at 
the Weyl nodes. The sign change in the Dirac mass sig- 
nals a quantum Hall transition. For small \k z \, where the 
mass is negative, the system is still in the topologically 
trivial regime. The 2-dimensional systems correspond- 
ing to momenta outside the Weyl nodes, however, are 
in topologically nontrivial quantum Hall state. Chiral 
surface modes appear on any surface that is not perpen- 
dicular to z for each value of k z between the Weyl nodes. 
This restriction gives rise to so-called Fermi arcs. 



A. Effect of pairing on nodes 

When superconductivity is turned on, the Hamiltonian 
decomposes into two copies of itself, acting on Bogoliubov 
quasiparticles rather than electrons. This is due to the 
fact that superconductivity splits each electronic state 
into a particle-holc symmetric and particlc-holc antisym- 
metric state with an energy separation ~ 2 |A|. For the 
Hamiltonian H— , these new states correspond to g?£ and 
/jj defined in Eq. (p"5|) . 

Remarkably, the system does not develop a supercon- 
ducting gap, but rather each Weyl node splits into two 
separated Bogoliubov Weyl nodes of opposite particle- 
hole symmetry (see Fig. [2Jb)) . Both Bogoliubov Weyl 
nodes have the same chirality, which is inherited from 
the initial electronic Weyl node, and in this sense half 
of the topological charge of the initial Weyl node. The 
particle-hole symmetric and particle-hole antisymmetric 
subspaces are decoupled. For fixed k z , each subspace de- 



scribes a spinless p + ip-supcrconductor, which is known 
to have both a topologically trivial ( "strong pairing" ) and 
non-trivial ( "weak pairing" ) phased The transition be- 
tween the strong and weak pairing phases is still marked 
by the Weyl nodes that separate trivial state at small \k z \ 
from a non-trivial state at large \k z \. 

If superconductivity is further increased, the distance 
between the particle-hole symmetric and antisymmetric 
Weyl nodes grows. This increases the topologically non- 
trivial momentum range for one of them (the 7£t A sub- 
space), and shrinks it for the other one (the %Z A sub- 
space). For HZ A , the Weyl nodes are pushed back to- 
wards k = (0, 0, n/d) T , where they annihilate. The cor- 
responding subspace is left in the topologically trivial 
insulating state. The Weyl nodes of H.Z. A move towards 
k = 0. After their annihilation at the origin, they leave 
the entire Brillouin zone of this subspace in a topologi- 
cally nontrivial insulating state. The motion of the Weyl 
nodes upon increase of |A| is followed in Fig. [2j Specifi- 
cally, subfigure (b) depicts the situation where both sub- 
spaces have topologically trivial and nontrivial momenta. 
Subfigure (c) corresponds to an order parameter ampli- 
tude |A| large enough to trivially gap out the HZ A sub- 
space, while the "Hi A subspace still has Weyl nodes. Sub- 
figure (d) corresponds to even larger |A|, such that all 
Weyl nodes have annihilated. 



B. Majorana surface states 

For the topologically non-trivial momentum range of 
k z , surface states are expected. We model a surface per- 
pendicular to y by replacing m and |A| by some smooth 
functions of y with m(y), |A|(y) = const, for y < 0, and 
m(y), |A|(y) — s- for y —> +oo (which realizes a trivial 
insulator equivalent to the vacuum). The Hamiltonian 

nt A = Vf (k x a y + i-^aA + M+ A {k z ,y) a z (23) 

indeed has eigenstates 

=± r ei*«^'»M» ( e ;:/ 4 4 

(24) 

which arc normalizable and exponentially localized at the 
surface only for momenta k z with M+ A (fc z ) > inside 
the sample, as anticipated. Af is the corresponding nor- 
malization factor. The dispersion of the surface state is 
linear, E = vpk x /2. 

The "Majorana-ness" of this state can be understood 
by counting. In particular, recall that the Nambu con- 
struction, Eq. nominally doubles the number of 
components of the fermionic fields. This implies that 
ipk and ip-k are not independent. Thus corresponding to 
Eq. (|2~4"|) there is one state, i.e. one canonical (complex) 
fcrmion, for each momentum satisfying the localization 
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condition (k z ) > with, say, k x > 0, where the last 
condition is made to keep the states independent. Equiv- 
alently, we can divide this complex fermion into two real 
ones, and associate one real Majorana fermion with each 
k x , with no restrictions on k x . 

In conclusion, we find that the Bogoliubov Hamilto- 
nians "H^ A and HZ A essentially each describe half of a 
normal Weyl semimetal. They are subject to respective 
effective Zeeman gaps m± | A | . A pair of Bogoliubov Weyl 
nodes of opposite chirality exists if m c \ < m±|A| < m C 2. 
In the sense that two Bogoliubov Weyl nodes arise from 
a single normal Weyl one, each of the former carries 
half of the topological charge of the latter. However, 
this notion is tied up with the non-independence of Bo- 
goliubov states. The corresponding quasiparticles are 
characterized by their chirality and particle-hole symme- 
try rather than just chirality in the non-superconducting 
case. The Bogoliubov Weyl nodes are located at mo- 
menta k = (0, 0, k z ) T , with 

A/ ±A (fc z ) = . (25) 

The vanishing of M±(k z ) reflects a topological transi- 
tion, regarding the quasiparticles as two dimensional ones 
parametrized by k z . For any k z that satisfies 

M± A (k z ) > , (26) 

the respective Hamiltonian T-L^ A maps to a topologically 
non-trivial spinless p+ip superconductor and has a chiral 
surface mode. The latter describes a Majorana particle 
with linear dispersion perpendicular to z and the surface, 
E = VFk ■ (ej_ x z) /2. The spin of the surface mode is 
locked to the direction of propagation. Negative values 
of M^ A (fc z ), on the contrary, correspond to a trivial in- 
sulator. 



C. Topological phase diagram 

The Hamiltonians can separately be tuned from 

a topologically trivial to topologically non-trivial state by 
changing both m and |A|. In our model, these two pa- 
rameters can be tuned separately, although a finite mag- 
netization can in principle affect the proximity induced 
superconductivity in the TI/SC interfaces. An analysis 
of the effect of the magnetization on the gap (and of the 
superconductivity on the magnetic ordering) requires a 
theory of the superconducting mechanism, like BCS the- 
ory, which goes well beyond the treatment here and we 
think deserves a separate study from this manuscript. 

If we exclude a substantial magnetic field due to mag- 
netic impurities, they could most importantly affect the 
superconductivity in the TI/SC interfaces by an exchange 
coupling of the nearby superconductor layers to the mag- 
netic impurities. A weakening of the superconductivity 
there would in turn diminish the proximity effect at the 
interfaces. Since however the superconductor layers are 




FIG. 3: Phase diagram of a Weyl superconductor in a TI/SC 
heterostructure as a function of Zeeman gap m and proxim- 
ity induced superconducting order parameter amplitude |A|. 
The values of m c i and m C 2 are set by the tunneling ampli- 
tudes between the surface Dirac layers of the heterostructure 
depicted in Fig.[T] see Eq. (|10p . Each phase is characterized by 
rib, the number of pairs of bulk Bogoliubov Weyl nodes, and 
n s , the number of two-dimensional Majorana surface modes. 
The phases are labeled as nt,n s . The black dots locate the 
different subfigures of Fig. [2] 



separated from the magnetic impurities by the interfaces, 
the exchange coupling of the SC to the magnetic impu- 
rities is certainly much weaker than the Zeeman term in 
the interfaces. It is therefore reasonable to neglect the 
explicit dependence of A from m. Similarly, a proximity 
effect of the SC on the impurities potentially weakening 
their ferromagnetic ordering is disregarded. This approx- 
imation would break down if either the magnetic field of 
the impurities was important, or if the tunnel coupling 
between the superconductor and the magnetic impurities 
was strong. 

Neglecting their weak interdependence, a simultane- 
ous modification of m and |A allows to access a number 
of different phases with ri\, = 0, 1, 2 pairs of Bogoliubov 
Weyl nodes in the bulk and n s = 0,1,2 two-dimensional 
Majorana surface modes (potentially living in a restricted 
/c z -range). Fig. [3] shows the phase diagram as a function 
of m and |A|, the phases are labeled according to their 
values of (nt,n s ). The phase diagram is mirror symmet- 
ric for negative values of m. The tunability of Weyl su- 
perconductors (which can be mapped to four copies of a 
spinless p+ip superconductor per value of fc z , only two of 
which arc potentially topologically non-trivial) is similar 
to the one of superconducting quantum anomalous Hall 
insulators (which can be mapped to two topologically 
potentially non-trivial copies of a spinless p + ip super- 
conductor) , although the increased complexity leads to a 
richer phase diagram for Weyl superconductors.— Of par- 
ticular interest are the phases (0,1), which corresponds 
to a truly topological superconductor (in class D and thus 
with two-dimensional topological invariants similarly to 



a three-dimensional quantum anomalous Hall state), and 
(1,1), which is precisely half of a normal Weyl semimetal. 
We will analyze them further in the next section.— 

IV. VORTICES IN WEYL 
SUPERCONDUCTORS 

One of the most interesting features of two dimensional 
topological superconductors is that they may host a zero 
energy Majorana mode localized around a vortex. Collec- 
tions of such Majorana bound states allow for non-local 
storage of quantum information, which may make the 
stored information less sensitive to decoherence.— This 
motivates the analysis of vortices in Weyl superconduc- 
tors. We specifically discuss the behavior of vortices in 
the simplest (0, 1) and (1,1) phases of Fig. [3l which mini- 
mize the number of surface Majorana modes. The vortex 
physics in other phases is qualitatively similar, but may 
involve more Majorana modes. 

The suppression of superconductivity inside a vortex 
puts its core in either the (0, 2), (2, 2) or (0, 0) phase. For 
simplicity, we consider the core to be is in the trivially 
insulating (0, 0) phase. In general, the finite size of the 
core suggests that in any case the core cannot be sharply 
distinguished from a trivial state, in a full treatment. 
The (0, 0) state can always be realized for an appropriate 
choice of m and |A|. Nevertheless, our results are not 
affected by this assumption. Different values of m and 
|A| will at most change the number and/or direction of 
propagation of the interface Majorana modes. 

Because the heterostructure has one special direction, 
namely the z-axis along which the different layers are 
stacked, vortices parallel and perpendicular to z have to 
be distinguished, as depicted in Fig. [4] The qualitative 
physics of vortices can be understood by analogy to 3 He- 
A, which also is a Weyl superconductor^^ As pointed 
out in Ref. l32l . the momentum range in which Majorana 
bound states at a vortex exist is proportional to e„ • z, 
where e v is the direction of the vortex. In particular, a 
vortex perpendicular to z has no bound states. In the fol- 
lowing, we will analyze the behavior of vortices in TI/SC 
heterostructures in more detail. 
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FIG. 4: The two classes of vortices in Weyl superconductors. 
Subfigure (a) sketches a vortex along the superlattice axis z, 
with bound states along a tube through the whole sample. 
Subfigure (b) depicts a vortex perpendicular to z. Whereas 
there are no states bound to the vortex, the surface states can 
be used for Majorana interferometry (thick line). 



phase has interface modes for any momentum k z . We re- 
strict the discussion to | A | < m when all relevant physics 
happens in but the results can easily be generalized. 

Exploiting the cylindrical symmetry with respect to 
the vortex axis, we model the latter by a radially depen- 
dent Zeeman gap m and superconducting order param- 
eter A(f) = |A(r)|e^W), where ip(<j)) = ip - ($/$o)</> 
is the phase of the order parameter. The latter is now 
twisted due to the presence of a magnetic field ((f) denotes 
the angular coordinate) . The twist is proportional to the 
flux $ trapped by the vortex, which itself is quantized 
in units of <I>o = hc/2e = ir/e, as usual for superconduc- 
tors. The phase tpo corresponds to the superconducting 
phase without magnetic field. The radius of the vortex 
is considered to be R, and m(r) and |A(r)| are smooth 
functions interpolating between fixed values m and |A| 
for r > R, and \m\ < m c i, |A| =0 inside the core of 
the vortex. The magnetic field is B = B z inside the 
vortex and for simplicity assumed to vanish everywhere 
else. This gives rise to a vector potential 



A. Vortex along the superlattice axis 



At first, we turn to a magnetic field B applied along 
z, the stacking axis of the heterostructure. For modest 
field strengths, only few vortices are present, and interac- 
tions between vortices can be neglected. This situation is 
sketched in Fig. IHa). By assumption, the vortex core is 
in a topologically trivial insulating state. The boundary 
of the vortex is thus equivalent to an interface between 
a Weyl superconductor and vacuum and has 1 Majorana 
edge mode. If the Weyl superconductor is in the (1, 1) 
phase, this mode has a restricted range of momenta k z (it 
lives "between the Bogoliubov Weyl nodes"), the (0,1) 



A(r) 
Mr) 



A(r)e<f, , 
Br 

-2^ R 



BR 2 
2r~~ 



(27) 



@(r - R) 



in e</, direction, 
coupling k — > k 



that is taken into account by 
— eA in the Hamiltonian ([TJ. 



canonical transformation a x —> —o v ,o v — > a 1 



minimal 
After a 
the role- 
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vant Hamiltonian H_ in Eq. (fTTj) becomes 



H A 



|A(r)| e-^W) cr* H_ 



(28a) 



(28b) 



H A = M-(k z ,r)a*+v F ^_^ l l (28c) 

/ -e-**\ fid . .. , 
+ VF [e^ ){rM +le 

For any given k z , this Hamiltonian may be interpreted as 
two copies of a spinlcss p+ip superconductor threaded by 
a magnetic flux. By analogy, the vortex binds one Majo- 
rana zero mode per topological value of k z and per topo- 
logical sub-sector if it traps an odd number of flux quanta, 
and no zero mode for an even number of trapped flux 
quanta=2, Assuming that there is only a single topologi- 
cally non-trivial subsector, one can thus define a unique 
zero energy Majorana mode bound to the vortex. See 
appendix IB"! for more details. 

Physically, the Majorana bound state can be under- 
stood in terms of an Aharonov-Bohm like phase, a Berry 
phase and a geometrical phase for the Majorana surface 
states. Consider the topologically equivalent situation of 
a Wcyl superconductor with a tube-like hole along the z 
axis. Without a magnetic field inside the hole, we know 
that chiral Majorana surface states exist when k z is cho- 
sen in the range where the two-dimensional superconduc- 
tor is in the topological phase. Since the spin is locked to 
the momentum, the surface states pick up a Berry phase 
of 7r upon encircling the hole once. This shifts the zero 
momentum mode away from zero energy and can be in- 
terpreted as effectively antiperiodic boundary condition 
on the geometrical phase in order to counterbalance the 
Berry phase. If now a unit flux is threaded through the 
tube-like hole, the surface states pick up an additional 
phase of w. The latter derives from the winding of the 
order parameter phase, and is similar to an Aharonov- 
Bohm effect. It compensates the Berry phase and thus 
allows for zero energy bound states. Similar effects have 
also been discussed for confined magnetic flux tubes im- 
posed in 3-dimensional strong topological insulators! 34 ' 35 
For momenta k z which are in the topologically trivial 
range, of course, no bound states exist both with and 
without magnetic flux. Because the magnetic field van- 
ishes outside the vortex, the topological character and 
especially the existence of surface states is unchanged 
there. We thus conclude that a vortex with an odd num- 
ber of flux quanta traps a Majorana zero mode for every 
topologically non-trivial value of k z . 

In a more realistic model, the Majorana bound states 
do not form totally flat bands as a function of k z . The 
presence of a zero energy Majorana mode for odd-integer 
fluxes is however partially robust. As an effective model 
at lowest energies, we consider the zero energy band of 



Majorana modes as a function of k z . After transforming 
to Wannier orbitals, we obtain a set of Majorana bound 
states at different heights z, as depicted in Fig.[4ja). This 
Hamiltonian can be interpreted as a 1-dimensional chain 
of decoupled sites. Next, we introduce a small hopping 
along the chain, thus allowing the Majoranas to move up 
and down the vortex tube. In dimensionlcss units, their 
dispersion is given by 



E 



cos(k z ) . 



(29) 



Therefore, zero energy Majorana bound states exist if 

1. a Majorana bound state can be defined for k z = 
±7r/2, i.e. M ±A (±7r/2) > 0, and if 

2. k z can take the values ±7r/2. 

The first condition is always fulfilled in the (0, 1) phase, 
but depends on the exact position of the Bogoliubov Weyl 
nodes in the (1, 1) phase. The second condition depends 
on the number of layers of the TI/SC heterostructre and 
the boundary conditions. Choosing for instance hard wall 
boundary conditions, one finds exactly one zero energy 
Majorana mode if the system has an odd number of su- 
perlattice layers, and no zero energy Majorana modes for 
an even number of layers. This result is quite natural, 
in a weak tunneling picture. Majorana states in a pair 
of layers can mix to form a Dirac fermion, moving away 
from zero energy. Only for an odd number of layers is an 
unpaired Majorana left behind at zero energy. 



B. Vortex perpendicular to the superlattice axis 

Viewing the vortex, as in the previous subsection, 
as a cylindrical hole enclosing a flux, the discussion in 
Scc. lIII Bl implics that a vortex perpendicular to the axis z 
of the heterostructure should also host Majorana modes. 
In our model, they run between the front and back sur- 
faces of the heterostructure on the side walls of the vor- 
tex. For a thin vortex, however, already a small coupling 
across the flux line is sufficient to hybridize and conse- 
quently gap out these two states. We thus recover the 
result of Ref. M- 

The hole comprising the vortex also introduces a new 
edge at the sample boundary (the circular ends of the 
cylindrical hole). The nearby surface states will rear- 
range in order to host the vortex and locally run along 
this new edge, as depicted in Fig. 0|b). While there are 
no states bound to the vortex, a special class of surface 
state paths allows for Majorana interferometry, depicted 
by thick lines in Fig. EJb). It may be interesting to study 
experimental measurements of interference in such struc- 
tures. 
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FIG. 5: Sketch of the anomalous thermal Hall effect in Weyl 
superconductors. The sample is shown from above. A thermal 
gradient VT is applied across the sample. The upper surface 
is at a temperature T> larger than the temperature T< of the 
lower surface. This leads to a net heat current from side A to 
side B perpendicular to the temperature gradient transported 
by the surface modes running around the sample. 



For concreteness, consider a temperature gradient VT 
imposed across the sample in the y direction. This leads 
to a net difference in the distribution of quasiparticles 
on the y = and y = L y surfaces. The result is an 
excess heat current Iq, in the x direction, which defines 
the thermal Hall conductance K xy , according to 



I Q = K xy |VT| 



(32) 



For small temperature differences between the surfaces, 
the excess heat current is obtain by differentiation, and 
we obtain 



1 k 2 B TT 2 T 

X" 2 3h ' 



dn F {v F k x ) 
dT 



(33) 
(34) 



V. EXPERIMENTAL PROBES 
A. Anomalous thermal Hall effet 

Because Majorana particles do not carry electric 
charge, a natural way to measure surface Majorana states 
is to detect their thermal transport. We again focus 
on a Weyl superconductor in the (0, 1) or (1, 1) phase. 
As there is only a single Majorana surface modes, these 
phases should exhibit only half of the thermal transport 
of a normal Weyl semimetal in the corresponding regime. 

On a surface perpendicular to y, thermodynamics can 
be calculated from the effective Majorana surface parti- 
tion function 



-S 



(30) 



S= *">n,fee,k* {-i^n + Vp k x ) ^ Un ,k x ,kz , (31) 

k x >0,k z 

where the operators ^ k k create excitations above 
the Bogoliubov vacuum. Note the restriction to k x > 0, 
which is because pairs of Majorana fcrmions at k x and 
— k x have been recombined into the canonical ^ fcrmion 
(c.f. Sec. Hill. 

If a thermal gradient VT is applied across the Weyl 
semimetal, each surface mode transports heat only in its 
direction of propagation. Therefore, the thermal gradient 
leads to a net heat transport perpendicular to VT, as de- 
picted in Fig. [5] This phenomenon is known as the ther- 
mal Hall effect. It has been proposed as an experimental 
signature of various other chiral edge states, for example 
in the spin Hall effect, the fractional quantum Hall effect 
or topological superconductors ! 28 ' 37 " — We presume any 
bulk transport to be parallel to the gradient, so that it 
can be separated from the surface contribution. In any 
case, the dependence upon field, density, etc. of any pos- 
sible bulk contribution would be very different from that 
of the surface one. 



with k z being summed over all topologically non-trivial 
values for the given phase of the Weyl superconductor 
(either (0,1) or (1,1)) and n F denoting the Fermi-Dirac 
distribution at the temperature T. Note that we have 
restored physical units such as Boltzmann's constant ks 
and Planck's constant h for concreteness. As expected, 
the surface of a Weyl superconductor has half of the ther- 
mal Hall conductance of a quantum Hall edge state per 
allowed momentum k z . This is not surprising because 
the thermal Hall coefficient is proportional to the cen- 
tral charge c of the surface modes, K xy = cK 2 k 2 B Tj (3h), 
similar to the heat capacity ! 28 ' 40 

Coming back to the Weyl superconductor in the (0, 1) 
or (1, 1) phase, the thermal Hall effect has an anomalous 
coefficient proportional to the distance 2k+ between the 
Weyl nodes defined in Eq. (fl"%|) . Concretely, the ther- 
mal Hall conductance is proportional to the length of the 
system in the z direction, K xy = n xy L Zl with 



i kWr a£ 



3/i 



(35) 



In the (0, 1) phase, where fc^ = n/d, each TI layer con- 
tributes the full Majorana quantum (1/2) ir 2 kgT/(3h) to 
the thermal Hall coefficient. Although thermal transport 
measurements are experimentally demanding, the higher 
dimensionality of the surface states in a Weyl supercon- 
ductor as compared to fractional or spin quantum Hall 
edge states hopefully tends to result in more approach- 
able experiments. 



B. Electrical transport 

As discussed in Sec. IIII1 the surface physics of a Weyl 
superconductor can be understood as layers of spinless 
p + ip superconductors stacked in momentum space along 
k z , with potentially associated edge states. In order to 
minimize bulk transport, we now specialize to the (0,1) 
phase. The surface of the Weyl superconductor is then 
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equivalent to just one non-trivial spinlcss p + ip super- 
conductor edge state per value of k z . In this phase, elec- 
tric transport experiments that have been proposed for 
p + ip superconductors can simply be transferred to Weyl 
superconductors. The general idea is to bring different 
samples with Majorana edge modes into contact. When- 
ever an interface has two edge modes running into the 
same direction, electrons can tunnel into the interface by 
decomposition in the two Majorana particles. These two 
Majorana particles can then be transported in parallel, 
giving rise to a one-directional electronic transport chan- 
nel along the interfaced In alternative setups, the two 
Majoranas can be separated and recombined with differ- 
ent Majorana modes, which leads to distinct signatures 
in conductance and noised The latter experiments are 
however less appropriate for Weyl superconductors where 
each surface has a large number of gencrically coupled 
Majorana modes at different values of k z . 



VI. SUMMARY AND CONCLUSIONS 

We have shown how a variety of gapless and/or 
topological superconducting phases can be achieved 
in superconducting-topological insulator superlattices. 
These phases are analogous to quasiparticle states of 3 He, 
The most interesting (0, 1) and (1,1) phases exhibit Ma- 
jorana surface states on some surfaces, and bound to the 
cores of vortices. Particularly in the gapless phases, such 
as (1,1), these Majorana states exist only for a range of 
momenta, k z , along the modulated direction of the su- 
perlattice. In such a case, no local (in z) description of 
the Majorana modes is possible, as opposed to the situ- 
ation in the (0, 1) phase, in which the Majorana modes 
can be modeled in terms of a real-space tight-binding 
Hamiltonian in the z direction, and the state can be con- 
sidered as a sort of stack of two-dimensional topological 
superconductors . 

It is hoped that the proposed structures might be ex- 
plored experimentally in the future. While we do not 
discuss materials in any detail here, we note that recent 
studies have shown that Cu^P^Se;?^— becomes a su- 
perconductor with x »s 0.14, while it is a topological 
insulator for x = 0, so that a super lattice with mod- 
ulated x might be a candidate realization of this pro- 
posal. Some ab initio modeling of such a superlattice 
would probably be useful prior to any experimental at- 
tempts. Alternatively, spin-triplet superconductors have 
recently been identified to exhibit Weyl superconducting 
phases as welL 43 ' 44 

There is significant scope for further theoretical study 
of Weyl and topological superconductors in three di- 
mensions. The Majorana surface states of (0, 1) and 
(1,1) phases are rather analogous to the "chiral surface 
sheaths" which occur in three-dimensional quantum Hall 
systems ; 45 ' 46 where interesting vertical transport, quan- 
tum interference, and universal conductance fluctuations 
have been studied, and it would be interesting to see 



how such phenomena translate to the superconducting 
case. We have also not touched on the Adler-Bell-Jackiw 
anomaly associated with Weyl points. This has been 
discussed recently for normal Weyl semimetals, where 
it may lead to anomalous magnetotransportii It is not 
obvious what the consequences are for Weyl supercon- 
ductors. One might also consider Josephson effects for 
currents along the z axis. We leave these questions for 
future work. 
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Appendix A: Gaplessness of inversion symmetric 
Weyl superconductors 

In this appendix, we discuss the qualitative behavior 
of Weyl semimetals under superconducting proximity ef- 
fect. We assume that either time reversal or inversion 
symmetry is conserved, while the respective other sym- 
metry needs to be broken for the system to exhibit Weyl 
physics. As has been mentioned in the introduction, the 
fate of the Weyl superconductor depends on which of the 
symmetries is conserved. For inversion symmetric, time 
reversal symmetry broken Weyl semimetals, the presence 
of a Weyl node at ko implies the presence of a Weyl node 
of opposite chirality at — fco- In Weyl semimetals with 
broken inversion symmetry, however, time reversal sym- 
metry guarantees the presence of Weyl nodes of equal 
chirality at ±/co- When a superconducting proximity ef- 
fect is turned on, the low energy modes at these two 
Weyl nodes mix. However, the superconducting correla- 
tions more precisely couple electrons on one Weyl node 
to holes on the other, instead of electrons to electrons as 
a more standard perturbation would do. This effectively 
inverts the chirality of one of the Weyl nodes. Conse- 
quently, the proximity effect in inversion symmetric sys- 
tems effectively mixes Weyl nodes of the same chirality, 
and no gap opens. In time reversal symmetric systems, 
on the contrary, the mixed Weyl nodes effectively have 
opposite chiralities, and a gap is to be expected. 

To be more concrete, we consider the effective low en- 
ergy theory of a Weyl semimetal, which corresponds to 
electrons living close to Weyl nodes. Since standard su- 
perconductivity couples electrons at momenta ±k, we fo- 
cus our effective model on two of the Weyl nodes located 
at momenta ±fco- For an inversion symmetric system, 
where the nodes are of opposite chirality, this already 
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describes a complete minimal model. For a time rever- 
sal symmetric system, the nodes at ±fco have the same 
chirality, and there must exist at least two additional 
nodes, say at ±fci, of the respective opposite chirality. 
Since these two pairs of nodes are decoupled, we can un- 
derstand the system as two copies of the following Hamil- 
tonian (|A4|) describing only two Weyl nodes. The latter 
thus allows us to decide on the presence or absence of a 
gap. 

As advertized, our result will only depend on whether 
the two initial Weyl nodes have the same or opposite 
chirality. We therefore assume one Weyl node to have 
positive chirality, Hi ~ 3 ■ k, while the second node is 
so far keep in a general notation, H2 ~ ± ■ k in order 
to tackle both the time reversal symmetric and inver- 
sion symmetric cases simultaneously. The electrons close 
to these nodes are described by the operators c 1 2 and 
c 2fccr> respectively. We furthermore measure the mo- 
menta relative to the respective Weyl nodes, such that 
the non-superconducting Hamiltonian reads 



k V t,k,u 

fc 

^ fc,0 



(Al) 



C 2,fc,t 

2,fc,4.y 



vp a ■ k 



±v F a-kJ fc.o 



can be written as 



fc.A fc.A 



w = 



vp a ■ k 



(A4) 



a(fc)* 1 CT + /3(fc)* • a 



a(fc) 1 CT + fi(k) -3 
^fv F a ■ k 



(A5) 



For an inversion symmetric system, where the lower 
(plus) sign applies, the diagonal is proportional to the 
unit matrix in Nambu space. Superconductivity can 
therefore never open up a gap, but only shift the Bo- 
goliobov Weyl nodes in momentum space. For time re- 
versal symmetric systems, the Hamiltonian is however 
generically gapped. As an example, we consider s-wave 
superconductivity. The latter corresponds to 



i,fc 



c t - , + h.c. 



t 2,-fe4 



(A6) 



^9 V x.fc 



c 1 - — c 1 - - ) + h.c. 

,t 2,-fe,4. 2,fe4 l,-fc,ly 



where we neglect the superconducting phase for simplic- 
ity, i.e. A = |A|. The total Hamiltonian then reads 

H = Y*t ( VF /-* . (A7) 

k 

The eigenvalues of this Hamiltonian are easily found as 



with %o = ( c i,fe,t' c i,k,V c 2,k,v c 2,k,i) ■ We assume that 
the superconducting part of the Hamiltonian only con- 
tains terms of the from 



H sc - h(k) c\ % 



a 2, — k,<j' 



h.c. 



(A2) 



which in particular includes s-wave and p-wave pairing. 
It is now useful to rewrite the non-superconducting part 
of the Hamiltonian as 



Hn 



vp a ■ k 






Tvf3 



fc,A 



where A 



(A3) 



We note 



(c - ,c - ,c' 

V l,fc,t l,fc,4- 2,-fc,i' 2,-fc,t' 

that the sign of the second Weyl node in the Hamilto- 
nian has to be reversed due to the inversion of creation 
and annihilation operators. This precisely corresponds to 
the effective inversion of the chirality of the second node 
under proximity effect, see above. Including a general su- 
perconducting term, the full Hamiltonian H = Hq + Hsc 



E = ±\ j{vp(T-k) 2 + 



IAP 



(A8) 



and the system is gapped as expected. Finally, consider 
adiabatically restoring inversion symmetry. Throughout 
this process, superconductivity ensures the system to be 
gapped. A time reversal symmetric, inversion symmetry 
broken Weyl superconductor can thus adiabatically be 
connected to the trivial state respecting both symmetries, 
and is therefore a trivial insulator itself. 



Appendix B: Expression of the zero energy 
Majorana bound state 

A vortex in a Weyl superconductor traps a unique zero 
energy bound state if it contains an odd number of flux 
quanta. To explicitly show this, let us first discuss the 
bound state for a simple limiting case where the alge- 
bra can be done explicitly, and then turn to the general 
solution. 

The limiting case is defined as follows. We assume 
that the Zeeman mass m is constant in the entire Weyl 
superconductor (and in particular takes the same value 
inside and outside the vortex). Moreover, we assume 
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that m Cl < m < m C2 , such that there is one momentum 
k z = k® with M_(r, kP z ) = everywhere. As follows from 
Fig. [3j we are then able to find a | A| = A outside the 
vortex such that only one subsector is topologically non- 
trivial, and expect a single zero energy Majorana mode 
bound to the vortex for this choice of |A|. In addition, we 
assume that there is only a single flux quantum inside the 
vortex. Outside the vortex, the Hamiltonian (|28[) reads 
for fc, = fc? 



= / d 2 r^o(f)W_V fe o(r) , 

lr>R 



H- 



Hi 



|A(r)|e^ «r 



Hb = vf 



|A(r)|e-^W)cr 2 H-b 
-ie-^\ d 



s*+ dr 



(Bla) 
(Bib) 
(Blc) 



-e"^\ /l 3 Si? 2 
+ V F i0 ~ II -— + ie 



e I<?> 



2r 



The order parameter amplitude |A(r)| goes to zero in 
the vortex core and takes the value |A(r)| = Ao far 
away from the vortex. For k z = k®, this Hamiltonian 
has two linearly independent normalizable zero energy 
bound state solutions, 



pouter J_ J_ - Sn dr' \A(k°,r')\/v P 
1 \ fll /~Z 



A/"' 



pouter _ J_ J_ - dr' |A(fe;,r')|/vj5 
2 ~~ A"" v 7 ?" 



Inside the vortex, i.e. for r < R, where 




(B2) 



(B3) 



Hb = v F 



only the state 



-ie-^\ d 



-ie l,p J dr 
-e-^\ fld_ . Br 
){rd^ +le ~ 



(B4) 



vrj inner 

2 ~ A"'" 



e -l/v F p R dr' (\A(k°y)\+eBr' /2) « ^gg-j 



is normalizable. The state that would be connected to 
is given by 



-l/v F f^clr' (\A(k°,r')\ + l/r' -eBr' /2) 








(B6) 

This state however diverges at the origin as •j/ , i nut ' r ~ 
1/r and is thus not normalizable. Consequently, there is 
only a single normalizable zero energy state bound to the 
vortex. Up to the normalization, it is given by 



* = ^ nncr Q{R - r) 



Q{r~R) 



(B7) 



When we consider a momentum k z close to fc° or change 
the Zeeman gap m a little bit, the system will stay in 
the same extended topological phase. There will thus 
always be a single zero energy bound state per topolog- 
ical momentum as long as there is no topological phase 
transition. We find that the Hamiltonian (|28p always ex- 
hibits two linearly independent zero energy bound states 
for r>R. For M+ A (k z ,r') > and MZ A (k z ,r') < at 
large r, they are given by 



pouter l _ p -S T R dr' M+ & {k z ,r')/v F 



V&°uter L P +fn dr ' MZ A (k z ,r')/v F 

* -A i— c 



i 
1 

\i e zrt > ) 

—i 
-1 
\ie i 4> J 



(B8) 



(B9) 



The bound state will be a superposition of these two 
states that connects to the normalizable solution for r < 
R. The special case considered previously corresponds to 



Ml* = -MZ A = \A(r)\ and f 2 ~ f+A - *z 
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